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ELECTRIC FIELD EFFECTS ON ION CURRENTS IN SATELLITE WAKES* ... 
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Small currents associated with satellite spin, dielectric conduction, or — 
trace concentrations of H + , can have a substantial effect on the potential 
of a satellite and the particle currents reaching its surface. The importance 
of such small currents at altitudes below about 300 km stems from the ex- 
tremely small 0 + currents impinging on the wake-side of the spacecraft.- 

The focus of the present study is the particle current on the downstream 
side of the AE-C satellite. Theoretical estimates based on a newly described 
constant of the motion of a particle indicate that accounting for small 
concentrations of H + remove a major discrepancy between calculated and 
measured currents. 


1. INTRODUCTION 

Many studies, both theoretical and experimental, have been made of the 
interaction between a satellite and the near earth plasma (refs. 1-11). The 
present study concerns charged particle current on the wake-side surface of a 
spacecraft in the earth^s ionosphere, where the vehicle is mesothermal; its 
speed V 0 exceeds the thermal velocity of plasma- ions but is much less than 
tne thermal velocity of the ambient electrons. 

In our analysis calculated currents are compared with those that have 
been observed by Samir et al. (ref. 10) on the Atmosphere Explorer C (AE-C). 
The AE-C experiments are weTT suited to our purpose, since its rate of spin, 
as well as the plasma densities, constituents and temperatures were known. 
Moreover, measurements were conducted at night, thus avoiding complications 
associated with active solar arrays. 

Measured satellite voltages were in the range V ~ 9-100, where the 
electron temperature 0 is in electron volts. These results exceeded 
theoretical estimates based on balance between ion and electron current on a 
conducting surface by a factor ~2.0-2.5. Theoretical surface potentials were 
also substantially less than expected from ion-0 + current balance at each 
point of a dielectric surface. In a previous work it was shown that ac- 
counting for either rotating currents of charge embedded in the dielectric or 
small concentrations of H + suffice to bring theory and experiment into con- 
formance (ref. 12). The calculated voltages reported in Reference 12 will be 
used in the calculations of current given in Section 3. 


* This work supported by Air Force Geophysics Laboratory, Hanscom Air Force 
Base, MA, under Co r ‘,,.< act No. F19628-82-C-0081. 
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2. ION CURRENTS 


^ The normal component of ion current density at a surface element loca t ed, 
at r s on the body is given by 
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v.n<0 


> » » > 
v.n f(r s ,v) d v 


( 1 ) 


where n is the outward normal at t’g. The distribution function f(r,v) at 
the phase space point r, v, satisfies, in general, the Vlasov-Poisson system 
of equations. For a perfectly absorbing body f satisfies the boundary 
condition- fir’s, v) = 0 for v.n > 0. 


Calculations of current to a satellite surface often use the assumption 
that ion currents to the satellite are given by the neutral approximation, 
which neglects the influence of electric fields. This assumption is quite 

good at the front (upstream) surface where ions reaching the satellite have- 

energies (~5 eV for 0 + ) substantially larger than electrical potential 
energies. The situation is less clear on the wake-side where electric fields 
may substantially enhance particle and current densities over the neutral 
approximation values. In the following paragraphs, this problem is. addressed 
by formulating bounds on j(?s) and and applying these bounds for an 

assumed, non-self-consistent model potential. 


The normal current density at a point r $ on the surface where the 
potential is V(r s ) can be written (m, = e = 1] 



Here v n = (v Q , fl 0 ) is the velocity of a particle at ? s , 

Po = .ft(r s ) (1 > v 0 > 0)> and % = < v > G) the velocity |t r - - on 

the trajectory that connects to the phase space point (v 0 , r s ). If 
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then 


j < j b * 2/f o (v) v 2 - V(r $ )j v dv d n , 

that is, jb is an upper bound on the normal ion current density at r s . 
Similarly the particle density satisfies 

1/2 » 
v dv d R 
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The bounds established herd require tbe inequality (3); that is, in terms of 
inside-out trajectories, neighboring orb-its emanating f rom r 0 with a given 
energy must diverge more at r = than at their point of origin. Although the 
inequality (3) appears to be a- reasonable assumption- for attractive poten- 
tials, the general conditions under which it applies have not been 
established-. 


for the following considerations, we take a Spherical satellite in the 
potential V(r,e) where r and e are spherical polar coordinates with polar axis 
in the direction of V 0 . We consider model potentials of the form 


V = -r~ 2 f(e) + V 0 (r) 


where V 0 (r) is a spherically symmetric potential. The asymptotic, far wake 
potential has in fact this form with (ref. 1) 


V 0 (r) * 0 


f(e) ~ a^ cos~^e exp(-M^ tan^e) 


where a is the radius of the satellite. For the discussion below the 
particular form of f(e) is arbitrary, however, and may be chosen to fit 
potentials near the satellite. 


The utility of the potential of the form in equation (6) is that a 
particle moving in it possesses three constants of motion. In addition to 
energy and the component of angular momentum about the polar axis the quantity 


C = -jj — f(e) , 


where L is the magnitude of the angular momentum about r = 0, is conserved 
along a part icl|, trajectory. This follows readily upon taking the scalar 
product of L = r x \Twith both sides of the torque equation (m^ = e = 1) 


= -r x VV 


taking account of equation (6) for V. There are fewer constants of motion 
than for V = V(r), since in the latter case the direction of L as well as its 
magnitude is constant. The dynamical constant C of equation (8) is a rigorous 
constant for potentials of the form (6), and should not be confused with the 

Whippirfref^lS^ Sam1 *' and Jew (ref * 14) and criticized by Laframboise and 


The effective potential for radial motion of a particle is now given by 


■eff< r > - V) * % 







With this V e ff we can solve for the orbits and evaluate the bounds on Ion 
current and particle densities at the satellite surface. The orbit equations 
are particularly tractable for the interesting case of particles which reach 
the poles-of the sphere, especially the wake-side- pole. For these particles 
the axial component of angular momentum is zero, the orbits are planar, and 
the solution of the equations of motion is reducible to quadratures. 

From 
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for the final direction e of a particle launched from r * a, e = n in the 
(Jirecti^n p 0 with speed v 0 . Here u 0 is the cosine of the angle between. 
v 0 and r 0 . 

Consider now a positive energy particle with initial coordinate 
Op » w. For attractive potentials, the orbits will appear as indicated in 
f vgure 1. 



# 0 = 7r 


Figure 1. Schematic diagram of inside-out orbits starting on the sphere r = 


The range of o integration for calculating the bounds of equations (4) and 
(5) is generated by variations of no between 0 and- 1. Here we have 
assumed that the particles- are not attracted -back to the surface- This affect 
however can only decrease the particle flux at the point of interest, and 

therefore does not influence the nature of. equations (4) and (5) as upper 

bounds. 

To proceed further, we now specialize to potentials of the form 


V 0 (r) . 0 
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which approximates the surface potentials given by Parks and Katz (ref. 12). 
The final angle at r » » on the inside-out trajectory corresponding to 
u 0 * 0 is given by 
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and the limiting angle for the outside-in orbit is 

a = TI - G 


(15) 


(16) 


When the velocity distribution remote from the spacecraft is Maxwellian, jj, 
and nj, are given by 


v 0 + v? 
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dv 
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where 
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3. RESULTS 

Numerical Integration- for the case M = Vq/vj - 8, |V(«)J = 15 vf 
gives for 0* ions 
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o' v T- 


Jb 
N.ev 


1.55. x 10 


,-10 


Hl 1 1 

~ = 2.64 X 10 _A1 
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The last ratio is to be compared with the density ratio 



= 6.22 x 10" 16 


for the case of zero electric field. Thus, although the effect of electric 
field may be to yield densities several orders of magnitude larger than those 
obtained in the neutral approximation, they remain extremely small in 
comparison with ambient ion densities, amounting to about 1 rrr 3 for 
Np ~ lO^ 1 nr 3 . By way of further contrast, the ion density ratio is 
also Small compared with the electron density ratio 


"e -15 -7 

7T - e. 5 = 3.06 X 10 ' . 


Let us now calculate the effect of electric fields on the density of 
hydrogen ion current striking the wake side of the satellite at e = it. The 
result of this calculation is intended to throw some light upon the issue 
raised by Samir and Fontheim (ref. 16) in their comparison between measured 
current ratios I (4) /I (90* ) and those calculated from Parker's model (refs. 
17,18). Parker's model is based on solutions of the Poisson equation in which 
ion densities are determined by particle tracking techniques. Only one 
species of ion is treated, however, and its mass is the mean mass of ions in 
the plasma. Samir and Fontheim contend that the two-to- three-order of 
magnitude discrepancy between measured and calculated current ratios might be 
removed (1) by properly treating the separate ionic components of the plasma, 
or (2) by considering the non-steady nature of the plasma environment. The 
results of our calculations, summarized in Table 1 below, indicate that proper 
treatment of the hydrogen component of the plasma suffices to remove the 
discrepancy between theory and experiment. 

Several observations are in order. First, Table 1 shows that both the 
measured ratio r and the estimated upper bound on the ratio exceed the ratio 
estimated from the neutral approximation, which one may reasonably expect to 
be a lower bound on r. Second, the measurement at 160 should give a somewhat 
greater ratio than would be observed at BO*; the neutral approximation 
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Table 1, Comparison of Measured and Calculated Current Ratios* 4 ) 
Case^ ,c ) 

n(h + )/n(oV‘ ,) 

CJ(180)/j(J(90)]J5 
Jb(180)/J(90) (f) 

[j(160)/J(90)]<g T< 
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(b) 

(cl 
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1.2fl(-3) 

l.Q?(-3) 

1.40(-3) 


8.8(-6) 

7.4{-6) 

?.8(-6) 

2.fl(-4) 

1.5{»2) 

i.3(-2) 

1.7U) 

0.48 

1.49(-2) 

2(-?) 

6.8(-3) 
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(e) 
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combined effect of electron and magnetic fields is not considered in this 
paper,. 


g|, . Fen-different ionospheric satellites in different env-ironments, for- 

example, in polar environments, the relative importance of the various 
physical effects may differ from that found, for AE-C in the conditions we 
investigated* Thus for satellites subjected to fluxes of energetic auroral 
electrons, field enhancement of wake-side collection could be a substantial 
effect. 
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